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ABSTRACT 

We present mid-infrared (MIR) observations, made with the TIMMI2 camera on the ESO 3.6 m telescope, 
toward 14 young massive star-forming regions. All regions were imaged in the band, and nine in the Q band, 
with an angular resolution of ~ 1". Typically, the regions exhibit a single or two compact sources (with sizes 
in the range 0.008-0.18 pc) plus extended diffuse emission. The Spitze r-Galactic Legacy Infrared Mid-Plane 
Survey Extraordinaire images of these regions show much more extended emission than that seen by TIMMI2, 
and this is attributed to polycyclic aromatic hydrocarbon (R4H) bands. For the MIR sources associated with 
radio continuum radiation (Paper I) there is a close morphological correspondence between the two emissions, 
suggesting that the ionized gas (radio source) and hot dust (MIR source) coexist inside the H II region. We 
found five MIR compact sources which are not associated with radio continuum emission, and are thus prime 
candidates for hosting young massive protostars. In particular, objects IRAS 14593-5852 II (only detected at 
17.7 //m) and 17008-4040 I are likely to be genuine 0-type protostellar objects. We also present TIMMI2 
A^-band spectra of 8 sources, all of which are dominated by a prominent silicate absorption feature (~ 9.7 /im). 
From these data we estimate column densities in the range (7-17) x 10^^ cm"^, in good agreement with those 
derived from the 1.2 mm data (Paper II). Seven sources show bright [Ne II] line emission, as expected from 
ionized gas regions. Only IRAS 12383-6128 shows detectable PAH emission at 8.6 and 1 1.3 fim. 
Subject headings: dust, extinction — infrared: ISM — stars: early-type — stars: formation 



1. INTRODUCTION 

Massive stars are known to be born deeply embedded 
within dense clumps of molecular gas and dust, with distinc- 
tive physical parameters: sizes of ~ 0.4 pc, masses of ~ 5 x 
IO^Mq, and dust temperatures of ~ 32 K (e.g.. lFaundez et akl 
120041) . The formation process of massive stars within these 
clumps is, however, still under discussion. Current theo- 
ries are coalescence (e.g., IBonnell et al.lI998t : lBonnell & Bate! 
12002). gravitational collap s e of single massive pre stellar cores 
(e.g.. iMcKee & Tanll2003t iKm mholz et all 120071) . and com- 
petitive accretion (e.g.. ' Bormell et al.n2004 ). At least up to 
masses corresponding to late O-type stars and early B-type 
stars, the observational evidence of accretion disks and molec- 
ul ar outflows favor s the collapse scenario (see the reviews 
bv lGarav & Lizanoi ri999: Beuther et al. 2007). The observa- 
tional study of the early stages of high-mass stars is difficult 
because of three main reasons: the rapid pre-main-sequence 
evolution of massive stars, confusion problems for individual 
studies since they form in clusters and are located at longer 
distances than low-mass stars, and the large columns of dust 
and gas (A^(H2) ~ 3 x lO^^ cm'^, iGaravet al.ll2007l) of the 
dense molecular clouds that harbor them, making them com- 
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pletely dark at visual wavelengths. Thus, the study of the en- 
vironment around recently formed massive stars is better per- 
formed through observations at infrared, millimeter, and radio 
wavelengths, where the extinction is much smaller 

To investigate the birth process and early evolution of mas- 
sive stars, we are carrying out a multiwavelength study of a 
sample of 18 luminous IRAS sources in the southern hemi- 
sphere thought to be representative of young massive star- 
forming regions. The goal is to understand the physical 
and chemical differences between different stages of evolu- 
tion. The objects were taken from the Galaxy- wide survey 
of CS(2^1) emission jBronfman et al.l 1996) toward IRAS 
source s with far-infrared colors typical of compact H II re- 
gions (IWood & Churchwellll989l) . We selected sources based 
primarily on the observed CS(2^1) line profiles; looking for 
self-absorb ed profiles cons istent with inward or outward mo- 
tions (e.g., Mardonesl[T998 '), and/or with extended line wings, 
possibly indicating the presence of bipolar outflows. In addi- 
tion, the sources were required to have IRAS 100 iim fluxes 
greater than lO"' Jy and to have declination 6 < -20°. The lu- 
minosity of the IRAS sources, computed using the IRAS fluxes 
and distances derived by L. Bronfman (2006, private commu- 
nication) are in the range 1 x 10'*-4 x 10^ Lq, implying that 
they contain at least one embedded massive star. The charac- 
teristics of the radio continuum emission at 1.4-8.6 GHz of 
the ionized gas and of the cold dust emission at 1 .2 mm associ- 
ated with the IRAS sources in our sample are reported, respec- 
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tively, in the first and second papers of the series (iGarav et alj 
120061 12OO7I hereafter Paper I and Paper II, respectively). 

Here, we report the results of mid-infrared (MIR) A^-band 
and 2-band imaging observations, and A^-band spectroscopy 
toward a selection of sources in our original sample. The 
MIR continuum emission was imaged within regions of typi- 
cally ~ 30" X 30" centered on the IRAS sources, with an an- 
gular resolution of ~ 1". In addition, Spitzer images at 3.6, 
4.5, 5.8, 8.0, and 24 /im were obtained from the public data 
of the galactic surveys Galactic Legacy Infrared Mid-Plane 
Survey Extraordinaire (GLIMPSE) and MIPS Galactic Plane 
Survey (MIPSGAL). The primary goals of this study are to 
compare the MIR observations wit h the rad io continuum im- 
ages of similar angular resolution (iPaper Ih . to determine the 
characteristics and physical conditions of the warm dust sur- 
rounding recently formed massive stars, and to investigate the 
stellar content within the massive and dense cores. 

In Section 2, we describe the MIR observations and data re- 
duction. Section 3 presents the images and spectra, the overall 
results and computed parameters of the sample. In Section 4, 
we present further discussion, and in Section 5 the conclu- 
sions of this study. 

2. OBSERVATIONS AND DATA REDUCTION 

The A^- and g-band images and A^-band spectra 
were obtained using the TIMMI2 mid-infrared camera 
jReimann et al.l 120001) mounted on the ESO 3.6 m telescope 
at La Silla, Chile. The camera uses a Raytheon 320 x 240 
Si:As impurity band conduction high background array. A 
pixel scale of 0.3" was used for A^-band imaging observa- 
tions and 0.2" for Q-band images. We performed long-slit 
low-resolution spectroscopy (70" north-south oriented slit, 
AA/A ~ 160) using a pixel scale of 0.45" and a slit width of 
1.2" or 3.0". 

The observations were carried out during 2003 May 23- 
25. The TIMMI2 filters used were Nl (Aeff = 8.7 ^m, AA = 
1.2 ^im), N11.9 (Aeff = 11.7 ^m, AA = 1.2 ^m), and Ql 
(Aeff = 17.7 fim, AA = 0.8 /im). We imaged 14 sources in at 
least one filter in the A^-band region; and nine of them in the 
Q-band region. A^-band spectroscopy (8-13 /im) was made 
toward eight sources. The standard nodding/chopping observ- 
ing technique was used to remove the strong and variable 
thermal background emission. For imaging, most sources 
were observed in the small-source mode (chopping in the 
north-south direction, nodding in the east-west direction) us- 
ing large throws (> 30"), in order to detect the extended MIR 
emission. For spectroscopy, in general, smaller chop throws 
were selected (~ 20"). MIR calibration standards stars — 
HD 108903,HD 123139andHD 169916— were observed in 
different modes (imaging and spectroscopic) every couple of 
hours for photometric flux conversion and spectroscopic cal- 
ibration. Table [T] summarizes all the observations: Columns 
1 and 2 indicate the IRAS source identification and galactic 
coordinates, respectively; Column 3 gives the TIMMI2 filter 
(designated by its effective wavelength Aeff or by the range 
8-13 /im for the case of spectroscopy); Column 4 indicates 
the final on-source integration time; Column 5 indicates the 
field of view (FOV; east-west amplitude x north-south ampli- 
tude for images, slit width x north-south amplitude for spec- 
tra); and finally Column 6 lists the noise level after calibration 
(images: noise level in 1 arcsec^, spectra: average error in the 
full FOV integrated spectrum). 

The images were reduced using our own IDL routines based 
on scripts originally developed by M. Marengo (2004, private 



TABLE 1 

Summary of the TIMMI2 Observations 
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8.7, 11.7 and 17.7 refer to NI, Nil. 9 and Ql images, respectively, whereas the 
8-13 /im range corresponds to long-slit low-resolution spectroscopy. 
^ FOV estimated depending on the observing mode. For the small-source imaging mode 
(nodding in the east-west direction, chopping in the north-south direction): FOV = nod 
throw X chop throw. For the classic imaging mode (both nodding and chopping in the 
noith-south direction): FOV = arr"ay .v-size (96" for and 64" for Q) X nod throw. For 
spectra: FOV = slit width X nod throw. 

^ For images, the value listed is the noise level in 1 arcsec" after flux calibration; for 
spectra, we give the average error in the spectrum that includes the whole emission. 

communication). Basically, the chopping and nodding pairs 
corresponding to a single chop-nod cycle were co-added in an 
individual frame, and then each useful portion in that image 
was bari-centered before the combination of all the frames. 
This prevented us from losing angular resolution due to little 
offsets between different chop-nod cycles. The final image 
for each source was produced co-adding all the observations 
weighted by the integration times. The flux calibration was 
made using the tabulated fluxes for the standard stars given 
on the instrument Web page (based on preliminary versions of 
the spectral models by Cohen et al. 1999). Finally, the images 
were Gaussian smoothed to a resolution of ~ 1.8", in order 
to enhance the sensitivity to diffuse extended emission and to 
match the angular resolution of the Sp/feer-GLIMPSE images 
(see below), allowing us to compare both sets of data in a 
consistent way. 

For spectra reduction we used the IDL scripts provided by 
Siebenmorgen et al. (2004), slightly modified in order to in- 
clude the extended emission in the spectra extraction . This 
method is based on the optimal extraction described by Hornd 
( 119861) . and consists of summing along the spatial dimension 
after weighting by a source profile, which is computed by 
collapsing the image along the dispersion direction. Wave- 
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TABLE 2 
Observed Parameters 
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Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes. and arcseconds. The position is from the brightest component when 
the source is multiple; astrometry was performed using the Spifzei'-\KAC images. fP"'"' is the flux density calculated on compact sources using 2.0" radius aperture photometiy. F'"' is the 
flux density calculated on a region enclosing the total emission from each source. Icr enors are given, taking into account statistical and calibration uncertainties. 



^ The angular size 9 is the deconvolved FWHM, computed as = \/ FWHM- — FWHMsiti - , where FWHM is determined on the A'-band image (1 1 .7 /^m when available, 8.7 /^m otherwise) 
for the compact source, and FWHMstd is the averaged value of the standard stars in the corresponding filter. 
^ Angular size computed at 17.7 /^m. 
^ 3(7 upper limit. 

^ From lBrooks et alj (20031) , assuming a flux unceitainty of 10%. 



length calibration was done mainly using the tables given in 
the Web page. The same procedure was applied to the target 
and the correspo nding standa rd star, whose theoretical spec- 
trum (models by lCohen et alJ [T999i) was used to calibrate the 
final target spectrum. 

We complemented our observations with MIR data ob - 
tained by the Spitzer Space Telescope d Werner et al.l l2004h . 
Images at 3.6, 4.5, 5.8, and 8.0 fim were downloaded 
from the public data of the GLIMPSE Legacy Program 
(iBenjamin et al.l 12003 ^. that surveyed the plane in the four 
bands of the Infrared Array Camera (IRAC; Fazio et al. 2004). 
We also obtained images at 24 /im from the recently deliv- 
ered m osaics of the MIPSGAL Legacy Program ( Carey et alj 
l2009h . carried out using the Multiband Imaging Photome- 
ter for Spitzer (MIPS; Rieke et all [20041) . The 18 sources 
of our multiwavelength survey are within the il,b) range 
covered by GLIMPSE/GLIMPSE II and MIPSGAL, except 
IRAS 12383-6128 which has only data at 4.5 and 8.0 /im. 
The angular resolutions of the used Spitzer filters, calculated 
from the point-spread function (PSF) archives given in the 
telescope Web page, are 1.81", 1.76", 1.84", 2.14", and 
5.81" at 3.6, 4.5, 5.8, 8.0, and 24 fim, respectively. 

3. RESULTS 

3.1. TIMMI2 Images 

We detected MIR emission toward the observed regions 
in all the observed filters. Figures [T] and |2] show gray-scale 
images of the smoothed emission in the V-band (8.7 /im or 
11.7 /im) toward 12 regions. We do not show images of the 
emission toward IRAS 14095-6102 and 16524-4300, which 
have a low signal-to-noise ratio (S/N). In seven regions, we 
detected a single MIR source. Most of these have a cometary- 
like morphology, exhibiting a bright compact component at 



the head and a weak trailing extended emission. In three 
regions (IRAS 15502-5302, 17016-4124, and 17008-4040) 
we detected the presence of two distinct MIR objects, and in 
two regions (IRAS 16128-5109 and 17271-3439) the MIR 
emission has a complex morphology. The morphologies of 
the sources detected in both the N and Q bands are very 
similar The only exception is in the IRAS 14593-5852 re- 
gion, where we detected a source at 17.7 /im (14593-5852 II) 
which is not seen in the N filter (see Figure[3]i. 

Among the compact components, no one is strictly point- 
like: whereas the average FWHM size of the standard stars 
is 0.9" in the 8.7 /im filter, 1.2" at 11.7 /im, and 1.3" at 
17.7 /im, the compact components have FWHM sizes > 1.4" 
and > 1.5" at 11.7 /im and 17.7 /im, respectively. However, 
we do not discard the presence of point sources in a few cases, 
due to the uncertainties in the FWHM determination, which 
was made using a radial profile fitting. 

According to Saviane & Doublier (20Q5|), the pointing ac- 
curacy of the ESO 3.6 m telescope is ~ 5" which prevents 
a direct comparison of the MIR and radio continuum spatial 
distributions. We used the Spitzer-IRAC data to determine 
accurate astrometry of our images. The bright compact part of 
the smoothed TIMMI2 image was fitted to the corresponding 
peak in the longest wavelength IRAC image available (usu- 
ally the 8.0 /im image). In most cases, the morphology of 
the emission exhibited in the IRAC and our images are very 
similar, allowing us to match the positions without any con- 
fusion. The offsets between the original TIMMI2 coordinates 
and the new adjusted positions are typically o:^ 11". Since 
the G LIMPSE point-source accuracy is ~ 0.3" jMeade et alj 
l2005h . we estimate that our final TIMMI2 coordinates have 
an error < 1", produced by centering uncertainties (e.g., cor- 
responding regions not compact enough or peaks not well de- 
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(a) IRAS 12383-6128 imag: 11.7)im 
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Fig. 1 . — TIMMI2 W-band gray-scale images, overlaid with radio continuum ATCA contours. Source names are showed in the upper left corner, whereas the 
upper right comer indicates the specific TIMMI2 filter used ( 1 1 .7 or 8.7 /im) and the frequency of the ra dio c ontours (4.8 or 8.6 GHz). The images are displayed 
in logarithmic scale, over the full flux range. Contour levels are basically the same than those used in lPaper 1 
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Fig. 2. — Same as in Figure [T] (b) For IRAS 17008^040, no radio continuum emission was detected toward objects I and II up to a level of 1.5 mjy at 4.8 
GHz. 
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fined). 

Figures [T] and [2] also sh ow contour maps of the radio con- 
tinuum emission dPaper J, 4.8 or 8.6 GHz) overlaid on the A^- 
band images. These figures show that in most sources the spa- 
tial distribution of the radio and MIR emissions are highly cor- 
related, both in extended (IRAS 12383-6128, 13291-6249, 
14593-5852, 15502-5302, 16128-5109, and 17009-4042) 
and compact structures (the same sources, together with 
IRAS 15520-5234, 16458-4512, and 17016-4124 I). This 
suggests an intimate association between the ionized gas and 
the warm dust within these regions (see further discussion 
in Section |4.2| i. Although in five (of 10) compact sources 
the corresponding radio and MIR peaks are slightly shifted, 
the displacement never reaches more than 2± 1" and therefore 
does not affect the global morphological match. Since these 
shifts are small, they probably can be explained uniquely by 
positional uncertainties of the data. The exceptional cases in 
which an MIR source is not associat ed w ith radio continuum 
emission will be discussed in Section l4~4l 

Table |2] lists the source properties computed directly from 
the TIMMI2 images (positions, angular sizes, and flux densi- 
ties). For each source, the given position corresponds to the 
centroid of the compact component; the brightest compact ob- 
ject is taken as a reference when the region is complex. The 

angular size was calculated as 6* = y/FV^HNf - FWHMstd^, 
where FWHM is the full width half-maximum of the com- 
pact component and FWHMstd is the averaged value of the 
observed standard stars, both calculated by fitting the radial 
profile on the A^-band image (1 1.7 /im when available, 8.7 /im 
otherwise). We computed two flux densities: a "point-source" 
flux density F''°'"^, calculated on compact sources using a 
2.0" radius aperture photometry (diameter of about twice the 
smoothed resolution of 1.8"), and a "total" flux density 
calculated on a region (circular in almost all cases) enclos- 
ing the whole emission from each source. Photometric error 
was computed adding by quadrature the statistical noise of 
the image and the flux calibration error, which is mainly de- 
termined by the atmospheric fluctuations. We estimated this 
contribution as the standard deviation of the calibration fac- 
tor along the three nights of observation. The resulting total 
photometric error is ~ 10% for the 8.7 and 11.7 /im filters 
and 22% -34% for the 17.7 /im filter, which is affected by 
the presence of many atmospheric features within the filter 
wavelength range. 

Within the 14 regions observed with TIMMI2, we can iden- 
tify at least 15 simple bright compact sources with A^-band 
deconvolved FWHM sizes in the range 0.7" — 3.3" {0 in Table 
|2]l. Using the kinematical distances, the physical diameters of 
these objects are in the range 0.008-0.18 pc, with an average 
value of 0.03 pc. The physical sizes of the compact MIR com- 
ponents associated with radio continuum sources are in good 
agreement with each other, when measured at a comparable 
angular resolution (at 8.6 GHz). This supports the idea that 
the radio and MIR emission trace similar structures. 

3.2. Spitzer Images 

Spitzer and Q-band images of all sources are shown in a 
dedicated Web page ^ . In this Web page we further discuss 
the characteristics of the MIR emission toward each of the in- 
dividual regions of the survey and summarize what is already 
known about them in the literature. In general, the IRAC im- 
ages show similar structures to those seen in the TIMMI2 im- 
ages, although the former are much more sensitive to diffuse 



TABLE 3 
Spitzer FLUXES 



Source 


(mJyT 


f lot 

(mJyT 


f tot 

(mJyT 


/rtot 

8.0 fj,\i\ 

(mJy) 


(mJy) 


12383-6128 


a 


460 


a 


7980 


a 


13291-6249 


1140 


> 1990'' 


6290 


> 8370'' 


> 32400" 


14095-6102 


158 


219 


1770 


4410 


69300'' 


14593-5852 1 


255 


411 


1670 


4660 


> 32100''= 


14593-5852 11 


< 33' 


<44' 


< 347'^ 


< 950' 




15394-5358 


878 


5050'' 


> 7710'' 


9840'' 


46900'' 


15502-5302 I 


1160 


> 2380'' 


11000 


> 13300'' 


> 27800'='= 


15502-5302 11 


111 


174 


802 


2090 




15520-5234 


315 


936 


2760 


6560 


556000'' 


15596-5301 


128 


352 


1040 


1870 


44900'' 


16128-5109 


2860 


3950 


17700 


53300 


> 12300'= 


16272-4837 


87 


317 


432 


389 


6040 


16458-4512 


188 


423 


1890 


4530 


76000'' 


16524-4300 


366 


955 


2940 


5410 


65800'' 


16547-4247 


119 


564 


988 


913 


49900'' 


17008-4040 1 


557 


9640'' 


23700'' 


45100'' 


> 49500"= 


17008-4040 11 


27 


53 


359 


883 




17009-4042 


470 


1560 


5920 


> 12800'' 


> 63800' 


17016-4124 1 


71 


442 


1010 


2040 


195000'''= 


17016-4124 11 


25 


48 


242 


677 




17158-3901 


488 


528 


3000 


8350 


> 54300' 


17271-3439 


1020 


1430 


8440 


24600 


> 5130' 



' |http : / /www ■ mpif r-bonn . mpg . de/staf f /emorales/mir_ 



Note. — f "" is the flux density calculated on the same region as for 
T1MM12 photometry. Uncertainties are estimated to be ~ 20% for the non- 
saturated IRAC bands. 
" Out of GLIMPSE/MIPSGAL coverage. 

'' Saturated; lower limit obtained integrating the emission in the TIMM12- 
defined region. 

' Saturated; lower limit obtained using 35" radius aperture photometry. 
'' Saturated; recovered flux using the imageworks PSF fitting tool (see the 
text). 

"= Double sources not resolved. Flux includes components 1 and II. 
' Upper limit obtained integrating the diffuse emission in the TIMMI2- 
defined region. 

extended emission, at 5.8 and 8.0 /im. To illustrate this, we 
show in Figures [3] and |4] two examples of three-color images 
made combining the individual 3.6, 4.5, and 8.0 p,m IRAC 
bands, together with the corresponding and g-band images 
from TIMMI2. These two regions are those containing the 
most promising hi gh-m ass protostellar object (HMPO) candi- 
dates (see Section l4!4b . The diffuse emission seen in the 5.8 
and 8.0 pm IRAC filters, not detected with TIMMI2, is not as- 
sociated with radio continuum emission from ionized gas and 
probably contains a considerable contribution from polycyclic 
aromatic hydrocarbo n (P AH) emission bands at 6.2, 7.7, and 
8.6 /im (see Section |43] |. In addition, due to its higher sen- 
sitivity and shorter wavelengths available, the IRAC images 
reveal that in most cases the TIMMI2 source is the brightest 
one of a cluster of MIR objects. 

Some of the three-color images show the presence of IR 
features that appear greenish (the so-called green fuzzies), 
most of which are extended and not associated with radio 
continuum emission. This type of emission is found toward 
IRAS 15394-5358, 15520-5234, 15596-5301, 16272-4837, 
16547-4247, 17008-4040 (Figure S), and 17016-4124. 
Beutheretal. (2005) and Rathborne et al. ( 2005) first sug- 
gested that the bright and extended IRAC 4.5 /im emission 
might trace shocked gas by outflowing material ramming in 
the ambient interstellar medium. Indeed, numerical simula- 
tions predict that the contribution of H2 shock-excited line 
emission in the 4.5 pm band is an order of magnitude brighter 

individual/mir_individual ■ html| 
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o 
o 
o 
2? 



IRAS 14593-5852 

SpUzer-IRAC (3 colors) 



(b) imag: 8.7(j,m 

^^^^P^ cont:8.6GHz; 




imag: 17.7 jim 
,\ cont: 8.6 GHz 



lO UO lO lU I'i lU , , 

Right Ascension (J2000) 

Fig. 3.— IRAS 14593-5852. (a) Spitzer-IRAC image made using 3.6 (blue), 4.5 (green), and 8.0 /xm (red) bands, (h) T1MM12 8.7 fim gray-scale image, 
overlaid with ATCA 8.6 GHz contours, (c) TIMMI2 17.7 /im gray-scale image, overlaid with ATCA 8.6 GHz contours. 





Right Ascension (J2000) 

Fig. 4. — IRAS 17008-4040. (a) Spitzer-IRAC image made using 3.6 (blue), 4.5 (green), and 8.0 fim (red) bands. Also shown are ATCA contours (in blue) at 
2.5 GHz. (h) TIMMI2 11. 7 fim gray -scale image, wi th the locations of maser sites indic ated by small sym bols: OH main-hne by a white square tCaswell 1998), 
6.7-GHz methanol by a red curie tWalsh et all 19981) . and water vapor by a black cross <Forster & Caswelllil989ll . (c) T1MMI2 17.7 fim gray-scale image. 



than in the other IRAC bands dSmith & Rosenll2005h . Ob- 
servational support for this inte rpretation has been provided 
by rece nt Spitzer surveys (e.g.. Noriega-Cr espo et al.ll2004t 
ISmith e t al. 2006), which show that molecular outflows are 
particularly strong in the 4.5 /im band. We also found that 
most of the sources which exhibit green fuzzies (all but IRAS 
15596-5301 and 17008-4040) are embedded in large-scale 
8.0 /im extinction features known as infrared dark clouds 
(IRDCs), w hich host the earl iest stages of high-mass star for- 
mation (see Rathborne et al. 2005, and references therein). 

All the sources of our multiwavelength survey covered 
by MIPSGAL present very bright emission at 24 /im, quite 
saturated at the core of the PSF (with the exception of 
IRAS 16272-4837); in particular, the TIMMI2 objects lie 



just within these saturated cores. In some regions, interest- 
ing large-scale structures could be seen at 24 /im, hinting the 
presence of warm dust probably excited by our target sources 
beyond the FOV studied with TIMMI2. We tried to overcome 
the saturation and rescue the lost MIPS fluxes using the pub- 
licly available imageworks code written by T. JaiTett, which 
recovers the saturated pixels by applying a PSF model that is 
fit to the nonsaturated wings of the source. However, since 
most of our sources are embedded, the PSF wings are often 
contaminated by diffuse emission and the recovered flux was 
overestimated by the code. We therefore only consider valid 
those fluxes which are consistent with our spectral energy dis- 

'' [http : / /spider . ipac ■ caltech . edu/staf f / jarrett/irac/ index- 1 . 1 



8 



MORALES ET AL. 



tribution (SED) models, as described in Section 14711 they are 
listed in the last column of Table[3] For the remaining sources, 
we give lower limit fluxes, obtained by integrating the emis- 
sion of the PSF wings wi thin an aperture of 35" an d applying 
the aperture correction of Engelbra cht et al] (120071) . 

Table [3] gives also the results of the photometry obtained 
from the IRAC images. In order to compare properly with 
TIMMI2 data, we integrated the emission across the same re- 
gions as for the TIMMI2 photometry, even though the emis- 
sion in the IRAC images is usually more extended. This re- 
sulted, however, in a good correlatio n be tween MSX 8.3 /^m 
and IRAC 8.0 /im fluxes (see Section l^Tl i. indicating that the 
regions selected for photometry were appropriate. For the 
sources not observed with TIMMI2, we chose similar aper- 
tures (~ 6" -9" radius) enclosing the brightest source(s). Ac- 
cording to iMeade et al.l ([2005), magnitude uncertainties are 
estimated to be less than 0.2 for the GLIMPSE Catalog, cor- 
responding to a flux error of < 18%. Considering that our 
sources are not pointlike, we estimated our photometric un- 
certainty as 20%. For the saturated sources, we proceeded 
as for MIPS photometry, applying the imageworks code and 
considering the recovered flux only if it is consistent with the 
SED model; otherwise, we give a lower Umit flux integrated 
within the TIMMI2-defined region. 

3.3. TIMMI2 Spectra 

Figure|5]presents the 8-13 fim spectra of the eight observed 
sources, extracted within small angular apertures of ~ 4" -8" 
around the central positions of the spectral images. The north- 
south oriented slit was centered in the brightest part of each 
source. The spectra integrated in the whole emission regions, 
not shown here, exhibit practically the same shape and fea- 
tures, but logically with a different flux scale. We truncated 
the spectra at 8.2 /im due to calibration uncertainties at the 
lower wavelength edge. Error bars are plotted as filled re- 
gions around the main curve; they only represent the statisti- 
cal noise derived from the original spectral image and do not 
include other causes of error 

All the spectra exhibit a broad — and in most cases promi- 
nent — silicate band seen in absorption, at ~ 9.7 /im. Seven 
sources (all but IRAS 17008-4040 I, Figure lUf)) present 
[Ne II] line emission at ~ 12.8 /im, which is spatially ex- 
tended. This is consistent with the association of these sources 
with radio continuum emission of ionized gas. Column 2 
of Table |4] gives the [Ne II] line flux measured on the inte- 
grated spectrum of each source; the statistical uncertainties 
derived directly from the spectrum are < 2%. Emission in 
other two atomic fine-structure lines, [Ar III] at 9.0 /im and 
[S IV] at 10.5 /im, was observed toward IRAS 16128-5109 
(Figure |3d), but this is more evident in the integrated spec- 
trum). The [Ar III] is also present in IRAS 13291-6249 (Fig- 
ure |5jb)). The PAH features at 11.3 /im and 8.6 /im were 
detected in only one region (IRAS 12383-6128, Figure|5la)). 
The faint emission feature at ~ 9.7 /im in IRAS 13291-6249, 
15520-5234, and 17016-4124 1 (FiguresHa), (c), and (h), re- 
spectively) could correspond to the H2 0-0 5(3) line, but it is 
located in a wavelength range highly contaminated by telluric 
ozone features, which make noisier this part of the spectrum 
(note that this is in general not properly represented by the 
statistical noise shown in the spectra). We therefore label this 
line with a question mark in Figure |5] If present, the H2 0- 
S(3) line could be tracing a photodissociation region or the 
existence of shock events, since it can be either radiatively or 
coUisionally excited (see, e.g.. , Morris et al 2004) . 



TABLE 4 

Parameters Derived from TIMMI2 Spectra 



Source 


Flux [Ne II] " 




A'(H2)cold 


^(H2)l.2 mm 




(10"'^ 


(K) 


(10^2 




12383-6128 


9.7 


369 


8.9 


5.6 


13291-6249 


16.4 


759 


11.6 


9.0 


15520-5234 


7.1 


355 


13.6 


24.8 


16128-5109 


83.2 


233 


6.6 


9.1 


16458-4512 


10.1 


405 


13.8 


10.2 


17008-4040 I 


13.9 


339 


17.1 


21.4 


17009-4042 


31.1 


458 


13.8 


32.2 


17016-4124 I 


2.5 


. . . d 


... i 


31.3 



* Uncertainties derived from the spectra are < 2% (only statistical error). 

Best fit parameters; see Equation (T). Uncertainties are < 1%, but we 
consider them underestimated (see the text). 

Computed from Equation (Sj- Uncertainties are ~ 20%, only considering 
the error from the 1.2 mm flux. 

Low-quality spectrum. 

We fitted the continuum part of the TIMMI2 spectrum 
i.e., removing first the line features) using a simple model 
Pascucci et al.ll2004h which assumes that the MIR continuum 
flux arises from a small region of hot dust, whose emission 
can be described by a gray body with temperature T^, column 
density A^(H2)hot, and solid angle fid, surrounded by a larger 
cloud of cold dust with column density A^(H2)coid- Under these 
assumptions, the flux density as a function of frequency is 
given by 

= l]dBi.(rd)(l -e-'^-'^'^^^'-Oe"""'^^"''""'/^ , (1) 

where 

a^=RdgtJ,mHKi, . (2) 

Rdg is the dust-to-gas mass ratio, /i is the mean molecular 
weight per hydrogen molecule, nin is the hydrogen mass, and 
Ki, is the absorption opacity per mass of dust, whose tabulated 
values were taken from the Weingartner & Draine (2001) dust 
model, with Ry = 5.5. .Chapman et al. (2009) found that this 
model is more consisten t with the MIR extinction law o f high- 
density regions than the lWeingartner & Draind ( 12001b model 
with Ry = 3.1, because it includes larger dust grains and there- 
fore accounts for possible grain growth. The fitting procedure 
uses equation ([T]) with il^, T^, and A^(H2)coid as free param- 
eters. We constrain the parameter A^(H2)hot by imposing the 
condition that 

MH2)hot=-l/^A^(H2)cold , (3) 
V ^ 'c 

where ft^ is the solid angle subtended by the cold component. 
This condition, which arises by assuming constant density and 
spherical symmetry for both the hot and cold components, is 
imposed in order to overcome the degeneracy problem of the 
parameters A ^(H2 )hQt and fid, as was done for the SEDs fitting 
(see Section l4TTT l. We use as angular size of the co ld cloud 
that determined from the 1.2 mm emission dPaperlll) . 

Columns 3 and 4 of Table |4] list the parameters and 
A^(H2)coid obtained from the best fit to the spectrum extracted 
from the central part of each source. The best fit for each 
source is shown in Fi gure [5] as a dashed line. We assumed 
/?dg = 1/125 Id raind l2003i) and /i = 2.8 (mean molecular 
weight, adopting 10% abundance of He with respect to H). 
The uncertainties derived from the covariance matrix com- 
puted in the least-square minimization are < 1%, but we sus- 
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pect that they are significantly underestimated, since the fit- 
ting only considers the statistical noise in the spectra, which 
is also low. 

To assess the robustness of the derived values of the param- 
eters we investigated the effects produced by changes in the 
assumptions. First, we used the Weingartner & Draine model 
for the opacities with Ry = 3.1. Since for the TIMMI2 spec- 
tral range, the Ry =5.5 and Ry = 3.1 models do not present 
significant differences, we found that the fitting procedure is 
not sensitive to such change, and the adjusted parameters only 
differ by ~ 5% in Ta and ~ 3% in A^(H2)coid- 

The average optical depths between 8 and 13 /im inferred 
from the hot-component column densities A^(H2)hot derived 
using the fiducial model, are in the range 0.01 < {Ti,)g,-ij, < 
0.1, suggesting that the MIR emission from hot dust is indeed 
optically thin in these regions. Similar va lues for the opti- 
cal depths at 11.7 /im were derived by D e Buizer et al.l (120051) 
from the MIR flux densities of their sample of massive star- 
forming regions, which are mostly < 0.3. Hence, we also 
investigated the changes produced by assuming that the MIR 
emission is optically thin , i.e., using 

where A = A?(H2)hotf^d- Note that here the parameters A^(H2)hot 
and cannot be fitted simultaneously, but this allows us 
to avoid the constant-density condition (eq. |[3l). Under 
this assumption we derived almost the same column densities 
A^(H2)coid and temperatures than those given by the fiducial 
model (differences < 1%). 

Since in all eight sources the MIR emission is located at the 
peak of the millimeter core, the fitted value of the cold dust 
cloud column density , A^(H2)coid, can be compared with the 
column density, A^(H2)i.2min, derived from the peak value of 
the 1.2 mm flux density, ^j'^jmm llPapsr III) . Assuming that the 
dust emission is optically thin at 1.2 mm, then 

^peak 

A^(H2)l.2mm = 77-^ Ht^^ ' 

mm 

where 57b is the beam (24") and 7^ is the cold dust tempera- 
ture derived in Section 14. II W e use ki 2mm = 1 cm^ g~' from 
lOssenkopf & Hennind([T994l) , since the effect of coagulation 
and presence of dirty ice mantles on dust grains in dense cold 
cores is Ukely to be importan t at this wavelength, where the 
I Weingartner & Draind (120011) Ry =5.5 model might not rep- 
resent it properly (in particular, it ignores ice mantles). The 
column densities calculated using this relation are given in 
Column 5 of TableH) The column densities derived from both 
methods are in good agreement, within a factor of ~ 2, partic- 
ularly considering the uncertainties inherent to both methods 
(e.g., in the opacities and dust-to-gas ratio), as well as the 
differences in angular resolution between the TIMMI2 and 
1.2 mm observations. 

4. DISCUSSION 

4. 1 . Spectral Energy Distributions 

We constructed SEDs for all the 18 sources of our sur- 
vey, cofl ecting the flux densities at 1.2 mm (SIMBA data, 
iPaper It) , at 12, 25, 60, and 100 ^J.m (IRAS data), and at 8.3, 
12.1, 14.7, and 21.3 ^m, computed from the images of the 
Midco urse Space Expe riment (MSX) Survey of the Galactic 
Plane dPrice et alj|2001h . We also included the IRAC (nonsat- 
urated) and TIMMI2 fluxes hsted in Tables |2]and[3l sum- 



ming the flux of the two components for double sources. In 
most cases, a simple model consisting of emission from three 
uniform gray bodies with different temperatures and sizes was 
able to fit the SED. Since the regions exhibit generally com- 
plex MIR structures and probably present temperature gra- 
dients, this is indeed a coarse simplification. However, this 
model allows us to determine the average dust temperatures 
representative of each wavelength range. 

The MSX fluxes were computed just over the corresponding 
compact sources in the MSX images, in order to reduce the 
effect of the different resolution (~ 18") with respect to the 
IRAC data (~ 1.8"). Despite that difference, the SEDs show 
(Figure |6]) that the MSX fluxes are consistent with the IRAC 
and TIMMI2 fluxes (see, e.g., the good correlation between 
MSX 8.3 fim fluxes and IRAC 8.0 /im fluxes), indicating that 
the major contribution to the "compact" MSX flux arises from 
the brightest source(s) within its beam. The IRAS data, how- 
ever, have a much coarser resolution, so they include emission 
from a more extended region. Consequently, we often ignored 
the IRAS 12 /im flux for the fitting, and at 1.2 mm we con- 
sidered the whole flux to make it representative of the same 
region that IRAS accounts. Conservative flux uncertainties of 
20% were taken for SIMBA, IRAS, MSX, and IRAC data. For 
TIMMI2 data, we used the computed errors. 

The model assumes the presence of three spherically sym- 
metric components: a cold (c) extended component, a warm 
(w) inner component, and a hot (h) compact component. We 
consider them as radiating gray bodies and include the absorp- 
tion of the radiation by the enclosing components. The total 
flux density is then given by, 

F,=F^ + f: + F^ , (6) 

where 

= nhB,iTi,)( 1 - e-'^""'^ )e-(A'c+A'w)<T„/2 

F: = f)wB.(r„)(l -e-''»"-)e-'''"-/2 (8) 

F>f]eB.(rc)(l-e-(''/''°^') , (9) 

where the superscripts c, w, and h denote quantities for the 
cold, warm, and hot components, respectively; is the solid 
angle; T the dust temperature, and the H2 column den- 
sity. The cross section per hydrogen molecule, ct^, is given 
by Equation dill , where Rdf, = 1/ 1 25 an d fj, = 2.8, and we 
adopted the Wei ngartner & Draind d2001h = 5.5 model for 
dust opacities (see Section [33T l. For the fit we considered as 
free parameters j/q, (3, ilh, ^w, 7h, and T„. The parameters of 
the cold component, T^, fig, and were taken as fixed, and 
equal to the values listed in iPaper H . Further, we adopted as 
upper limits for ilh and fl^, the values derived, respectively, 
from the A^-band observations (Table|2]l and MSX /s-band ob- 
servations. Because the column densities for the hot and warm 
components are small, the expression {l-e~^'^") acts as a lin- 
ear term in A^, producing degeneracy between the parameters 

and ft. To overcome this problem, we constrained the A^h 
and A^„ parameters assuming constant density for all the tem- 
perature components, which implies Ah,w = \/^h.w/^cNc for 
the hot and warm regions. Given the simplicity of the SED 
model and the differences in the data used, a more realistic 
density distribution is unnecessary. 

Table |5] lists the parameters of the hot and warm compo- 
nents obtained from the best fit to the SED using the model de- 
scribed above. The FWHM angular sizes 6^ and 9„ were com- 
puted from the fitted soHd angles assuming Gaussian emission 
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Fig. 5. — TIMMI2 8—13 fim spectra of the central parts of the sources. The statistical noise is shown as a filled region around the main curve. The spectral 
features found are labeled and marke d wi th vertical dashed lines. We also show here the best fit model of the silicate absorption band (dashed curve) according 
to the procedure explained in Section [33l 
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Fig. 6. — SEDs and their corresponding fits as described in Section ITTI Flux symbols are a filled circle for 1.2 mm S IMB A, stars for IRAS, empty circles for 
MSX, squares for TIMMI2, and triangles for Spfeer-IRAC/MIPS (filled if they are saturated fluxes satisfying Equation jlOt ). 
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TABLE 5 






Parameters Derived from SED 




Fitting 






Source 


7h 










(K) 


(K) 


(") 


(") 


12383-6128 


361 


71 


0.43 


12.3 


13291-6249 


398 


93 


0.34 


3.9 


14095-6102 


342 


65 


0.26 


7.2 


14593-5852 


328 


117 


0.41 


2.9 


15394-5358 


432 


99 


0.46 


2.8 


15502-5302 


263 


103 


1.78 


6.6 


15520-5234 


276 


112 


1.17 


2.7 


15596-5301 


420 


63 


0.16 


7.2 


16128-5109 


323 


123 


1.02 


5.8 


16272-4837 


591 


75 


0.09 


4.1 


16458-4512 


330 


73 


0.38 


6.7 


16524-4300 


425 


68 


0.27 


9.0 


16547-4247 


504 


139 


0.15 


1.5 


17008-4040 


340 


98 


1.12 


6.2 


17009-4042 


400" 


132 


0.57 


5.2 


17016-4124 


400" 


96 


0.33 


6.1 


17158-3901 


331 


62 


0.49 


12.0 


17271-3439 


332 


181 


1.18 


3.7 



Note. — Best fit parameters, see Equa- 
tions Q and {8). Solid angles (f2) were con- 
verted to FHWMs (6) assuming Gaussian 
distributions. Statistical en'ors are about 
5%, 15%, 12%, and 30% for Th, T„, ^h, and 
0„, respectively. 

" Unstable fitting; reached upper limit of 
400 K. 

distributions (il = 1.133 x 9^). Covariance-matrix errors are 
about 5%, 15%, 12%, and 30% for Th, T„, 6^, and 9^, respec- 
tively, but they could be underestimated. Most of the hot dust 
temperatures fall in the range ~ 250-450 K, whereas warm 
dust temperatures fall in the range ~ 60- 150 K. 

Figure|6]shows the SEDs and their corresponding fits for all 
sources. In most cases the fit is excellent. For the objects with 
the largest cold cloud column densities {N^ > 3 x lO^-' cm"^), 
which gives rise to a very deep silicate absorption feature at 
~ 10 ^lm (e.g., IRAS 17271-3439, FigurelSr)), the inclusion 
of absorption in the model in the MIR range became difficult. 
Because of this, particularly unstable fits were obtained for 
IRAS 17009-4042 and 17016-4124 (FigureslSo) and (p), re- 
spectively). Hence, for these objects we set an upper limit of 
400 K for Th in order to obtain reasonable output parameters. 
The SEDs of the sources IRAS 16272-4837 and 16547-4247 
(Figures|6lj) and (m)) exhibit a peculiar tendency in the IRAC 
fluxes, resulting in high hot-dust temperatures (591 K and 
504 K, respectively). In these sources, the IRAC 4.5 /im flux 
density is similar (only slightly lower) to the 5.8 /im flux, 
which is higher than the 8.0 /im flux. This trend is consis- 
tent with the high absorption and the green fuzzies found in 
these regions. 

For each source, the best-fit SED model was used as a va- 
lidity check for the saturated IRAC/MIPS fluxes, which were 
recovered by applying the imageworks code (see Section lTTt , 
but were not used by the SED fitting. We list in Table [3] and 
show in Figure|6]as filled triangles, only those saturated fluxes 
that satisfy the condition 

where F^^ is the recovered flux at the saturated wavelength 



Ao = 4.5, 5.8, 8.0, or 24 /im, F™' is the flux density given 
by the SED model, is the number of data points being fitted, 
} are the observed fluxes, and /S.F\ represents the corre- 
sponding error (it was taken as 30% for the saturated flux). 
In other words, we request that Fj|iJ' did not deviate from the 

SED model more than the x^-per-data point value of the fit. 
Using this criterion, it was found that the imageworks tool had 
about 50% of success for embedded sources (e.g., seven out 
of 16 sources with saturated MIPS data). 

4.2. Warm Dust Regions and Association with Ionized Gas 

The high correlation between the spatial distribution of the 
MIR and radio emissions may suggest that the MIR emis- 
sion is due to the ionized gas. However, extrapolating the ob- 
served radio continuum flux of the optically thin H II regions 
jPaper I ) to the MIR, assuming a i/"" ' frequency dependence, 
we find that the free-free emission at MIR wavelengths is neg- 
ligible in comparison to the observed MIR fluxes (less than 
~ 6% at 8.0, 11.7, and 17.7 /im). We conclude that the MIR 
emission is tracing warm dust, loc ated either with i n the H II 
region or in a thin shell around it. iKraemer et alj ( 120031) al- 
ready found well-correlated MIR and radio distributions, and 
discarded the possibility that the emitting dust is distributed 
in a thin shell, since in that case the MIR emission would 
show evidence of limb-brightened morphology, which was 
not observed. None of the regions of our sample exhibit limb- 
brightened MIR emission. 

We conclude that the MIR emission of the sources asso- 
ciated with radio-continuum emission arises from warm dust 
mixed with the ionized gas inside the H II region. Although 
dust grains and gas particles are likely to be dynamically cou- 
pled in H II regions, the heating of the dust grains is dom- 
inated by the UV radiation from the exciting source(s) and 
the ionized nebula, and not by the colUsion al exchange of 
energy with the gas (see, e.g.. 'Sp itzeiill978 ). This makes 
possible the coexistence of a plasma with a kinetic temper- 
ature of ~ 7000 K and solid particles with a temperature of 
2± 250-450 K (Section l4TTTi. These dust temperatur es are con- 
sistent with the models by iNatta & Panagial (Il976h for dusty 
H 11 regions. 

4.3. PAH Emission 

The A^-band TIMMI2 images (and probably the 2-band im- 
ages too) are sensitive mainly to the emission of warm dust 
intimately associated with the ionized gas, or close to the 
exciting sources when there is no radio emission. They do 
not trace the presence of PAH, as indicated by the 8-13 /im 
spectroscopy toward seven of the eight sources of our sample 
(all but IRAS 12383-6128), where the emission bands in this 
wavelength range (the more prominent at 8.6 and 11.3 /im) 
are not seen. 

Within the H II region or very close to a massive proto- 
star, hard UV radiation destroys these molecules. For exam- 
ple, iPovich et alJ (l2007h calculated the photodestruction rate 
of PAHs by extreme-UV (EUV) photons and found that this 
destruction mechanism accounts for the absence of PAHs in 
the M17 H II region. Because this region is one of the most 
energetic in the Galaxy, we expect that for the UC sources of 
our sample, the PAH destruction edge is located much closer 
to the exciting star(s). IRAS 12383-6128, which is the only 
one that presents PAH emission bands at 8.6 and 1 1 .3 /im, has 
probably a softer UV spectrum than the others. 

It is possible that a part of the outer diffuse emission ob- 
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served in the TIMMI2 images is due to PAHs, but since 
the slit was centered in the bright compact MIR sources, 
the PAH radiation was not strong enough to be detected. 
This explanation is supported by other MIR-spectroscopic 
observations found in the literature, with larger apertures, 
toward some of our sources, in which PAH emission is 
detected. In particular, the IRAS LSR spectrum (aperture 
of 6' X 5') of IRAS 13291-6249 dJourdain de Muizon et alj 
119901) and the ISO spectrum (aperture of 14" x 20" in 
the wavelength range of TIMMI2) of IRAS 16128-5109 
(jPeeters et al. 2002) reveal the presence of PAH emission 
bands within the wavelength range of the TIMMI2 spec- 
tra (at 8.6 and 11.3 ^m). For IRAS 16128-5109, the 
PAH bands at 3.3, 6.2, 7.7, and 12.7 ^m are also present 
dPeeters et alj|2002l) . PAH emission was also detected to- 
ward other sources of our sample which were not observed 
spectroscopically by us: IRAS 14593-5852 (IRAS LSR^ 
Zavagno et al. 1992). IRAS 15502-5302 (ISO, Peeters et_al] 
2002t). IRAS" 15596-5301, and IRAS 17271-3439 (IRAS 



LSR. lJourdain de Muizon et alJ[T990 ). 

The latter result strongly suggests that the extended diffuse 
emission seen in most regions by the Spitzer-IRAC images 
beyond the radio continuum contours include an important 
contribution of PAH emission bands. The [3.6], [5.8], and 
[8.0] IRAC filters cover the PAH features at 3.3, 6.2, and 
7.7 fj,m, respectively. The [8.0] filter also covers the PAH band 
at 8.6 fim. This extended diffuse emission is, however, more 
clearly seen in the [5.8] and [8.0] IR AC images than in the 
[3.6] images. A possible explanation (Povi ch et al.l 12007) is 
that the 6.2 /im feature might have a greater fractional contri- 
bution to the [5.8] filter than the 3.3 /im band to the [3.6] filter 
(and the same for the [8.0] and [3.6] filters), or simply that the 
3.3 /im band is intrinsically fainter As for the TIMMI2 im- 
ages, the Spitzer MIR emission at the location of the UC H II 
regions is more likely to correspond to warm dust continuum 
radiation with little or any contribution of PAH features. 



4.4. High-Mass Protostellar Objects 

The high angular resolution MIR observations presented 
here revealed the presence of bright compact MIR sources 
in the vicinity of, but not coincident with, UC H II regions, 
which are likely to signpost HMPOs. We identify a to- 
tal of five new HMPO candidates within the regions ob- 
served with TIMMI2: IRAS 14593-5852 II (only detected 
at 17.7 /im, see Figure O, 17008-4040 I and II (Figure Hi, 
17016-4124 II (Figure |2|d)), and the compact component of 
IRAS 17158-3901 (Figure|2e)). 

In order to estimate the spectral types of these HMPO can- 
didates, we derived MIR luminosities from the measured A^- 
and Q-ba nd flux es, following a method based on the work by 
iDe Buizer et al.l ( 2000). This technique provides good lower 
limi ts to the true bolometri c luminosities of the stellar sources 
(see lDe Buizer et aLll2002l) . We note that estimates made us- 
ing the IRAS fluxes are inappropriate since they include the 
contribution from the exciting star(s) of the nearby UC H II re- 
gions, as well as from other IR sources within the large IRAS 
beam. First, a dust color temperature is computed assuming 
that the MIR emission from warm dust is optically thin, as 
we suggested in Section [33] Our SED fits are also consis- 
tent with optically thin emission from the hot-dust component 
(using the derived A^h, the optical depths at 10 /im are in the 
range ~ 0.01 -0.3). We adopt a standard expression for the 



TABLE 6 

Characteristics of the New HMPO 
Candidates 



Source 






Spectral Type 




(K) 






14593-5852 IP 


142 


54500 


08.5 


17008-4040 I 


219 


40900 


09.5 


17008-4040 II 


199 


2490 


B2 


17016-4124 II 


170 


6120 


Bl 


17158-3901 


169 


542 


. . . b 



Note. — Pai'ameters T^„\a^ and Lmir were com- 
puted from the ob.served N- and Q-band fluxes. The 
spectral types were derived from Lmir (and thus are 
lower limits) and the cahbration of iPanagia (|1973I) . 
" Lower limit W-band flux used. 
Spectral type later than B3. 

flux density (analog to Equation (|4|): 

f,=Aa,B,(W)e-''^"-/' . 



(11) 



Here, the factor A = Nh^ is constant, cr^ is given by Equa- 
tion (|2]i, and Nc by Equation (|5]l, but truncated to the max- 
imum column density derived from TIMMI2 spectra if it is 
higher than that value (17.1 x 10^^ cm"^, see TablelH, in or- 
der to moderate the effect of dust absorption (see below). The 
color temperature Tcoioi is derived from the ratio of the inte- 
grated fluxes at the A^- and Q-band (lower limit A^-band flux 
for IRAS 14593-5852 II): 



„lu;/(kX2T„io,) _ 1 F / \, \ ^ K 



ghc/(kXtTc„io,) _ I 



(12) 



The MIR luminosity is then obtained integrating an expres- 
sion of the form /v = A(T^B^(Tco\or) from 1 to 1000 /im us- 
ing the kinematical distance and the factor A to scale the flux 
density to the observed ones. Columns 2 and 3 of Table |6] 
give the computed color temperatures Tco\or and MIR lumi- 
nosities Lmir, respectively. We find that the derived spectral 
types, using the calibration of Panagia (1973) and the MIR 
luminosities, are earlier than B3 for four of our HMPO candi- 
dates (Column 4 of Table|6]l. 

The spectral types deri ved by us are ear lier than the typical 
spectral types obtained bv De Buizer et al.l (1200 0) using a sim- 
ilar method. A substantial difference between our technique 
and theirs is that we have considered the absorption of a cold 
cloud to model the observed fluxes, and corrected by extinc- 
tion when computing the MIR luminosities. We found, how- 
ever, that the employed method is quite sensitive to assumed 
column density Nc of the absorbing cloud, and therefore, in 
order to obtain strict lower limits for the spectral types, we 
recomputed the MIR luminosities under the assumption that 
there is no extinction (Nc = 0). This results now in two sources 
eai-lier than B3, namely, BO for IRAS 14593-5852 II and B2 
for 17008-4040 1, which thus might be considered as genuine 
HMPOs. 

Figures[3]and|4]present the whole set of MIR images for the 
regions containing the two bona fide O-type high-mass pro- 
tostar candidates: IRAS 14593-5852 II and 17008-4040 I. 
The Spitzer-IRAC images toward IRAS 14593-5852 (Fig- 
ure Oa)) show only diffuse emission toward the HMPO, 
bright in red (8.0 ^m) and thus probably corresponding to 
the environment's PAH emission, which dominates the large- 
scale structure. IRAS 14593-5852 II is then Hkely to be a 
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deeply embedded object. IRAS 17008-4040 1 is bright in all 
available MIR images, even saturated in the IRAC filters (Fig- 
ure Ua)). The H II region is displaced by ~ 30" to the west 
of the TIMMI2 sources and associated with extended diffuse 
emission at 8.0 /im. In addition to its high MIR luminosity 
and the lack of a detectable H II region (as indicated by the 
lack of both radio continuum emission and of the [Ne II] line 
in its A^-band spectrum), there is additional evidence support- 
ing the nature of IRAS 17008-4040 1 as a high-mass pro- 
tostar: (i) the Spitzer-IRAC image shows the presence of ex- 
tended filamentary green fuzzies (see Section [J!2] i to its north- 
east, suggesting outflowing activity; (ii) it is located at the 
center of a massive and dense molecular core (peak of 1 .2 mm 
emission. lPaper Ij) . as predicted by all the massive star forma- 
tion theories; and (iii) it is associated with different types of 
maser emission (Figure |4|b)), all of which are known to be 
associated with high-mass star-forming regions: OH masers 
at 16 65/1667 MHz (e.g., Forster & Caswe fl 1989; Caswell 
119981) and 6.7 GHz methanol masers (e.g.,|Walsh et al. 1998) 
just coincident wifli flie position of IRAS 17008-4040 1, and 
water vap or maser emission at 22 GHz located ~ 4" to its 
southeast ( Forster & Caswel]| |1989^. In particular, 6.7 GHz 
methanol maser s are only detected toward high-mass star- 
forming regions jMinier e t al. 2003). 

Finally, although the MIR luminosities given in Table |6] 
should only be considered as rough estimates, it is inter- 
esting to know where the HMPO candidates are located in 
the Lboi - Me,„ (bolometric luminosity - envelope mass) di- 
agram presented in the work of Molinari et al.l (2008, their 
Figure 9). They successfully used it to follow the pre-main- 
sequence evolution of massive young stellar objects, by clas- 
sifying a numerous sample of objects according to SED mod- 
eling and then co mparing with the pre dictions of the tur- 
bulent core model ( lMcKee&Tanll200l . We can do such 
exercise only for IRAS 17008^040 I, since it is the only 
HMPO candidate that is far enough from nearby UC H II re- 
gions and thus the mass of the associated 1.2 mm core can 
be considered as its envelope mass (neglecting also the sec- 
ondary object IRAS 17008-4040 II, which is much less lu- 
minous). Then, assuming Lboi = 40900 Lq (MIR luminosity) 
and Menv = 1200 Mq (Paper II), IRAS 17008-4040 I Hes in 
a zone dominated by MM-P o bjects, below that occu pied by 
IR-P objects in the notation of iMoUnari et alj (120081) . They 
suggested that the IR-P stage may correspond to the arrival 
of the protostar on the zero-age main sequence (ZAMS), and 
the MM-P objects represent an earlier stage. The location of 
IRAS 17008-4040 1 in the Lboi -A^env plot would then catego- 
rize it as a pre-ZAMS object. 

5. SUMMARY AND CONCLUSIONS 

We present MIR A^- and Q-b^nd imaging, made with the 
TIMMI2 camera, toward 14 luminous IRAS point sources 



with colors of UC H II regions, thought to be massive star- 
forming regions in early stages of evolution. We also present 
A^-band spectroscopic observations toward eight of these re- 
gions. Images from the Spitzer legacy programs GLIMPSE 
and MIPSGAL for the 1 8 sources of our multiwavelength sur- 
vey complement our data set. The main results and conclu- 
sions presented in this paper are summarized as follows. 

The morphology of the TIMMI2 emission toward most re- 
gions (10) is simple, exhibiting either one (seven cases) or 
two (three cases) compact components each surrounded by 
extended emission. The compact components have physical 
sizes in the range 0.008-0.18 pc. The S/^/fzer-GLIMPSE im- 
ages of these regions show even more extended diffuse emis- 
sion than that seen by TIMMI2, most likely due to emission 
from PAH bands. 

We find that the MIR emission traced by TIMMI2 corre- 
sponds to optically thin radiation from hot dust. The hot- 
dust temperatures, derived from fitting the SEDs, are typically 
~ 250- 450 K. 

We find that the spatial distribution of the MIR emission 
and radio continuum emission (Paper I) are highly correlated, 
suggesting that the ionized gas (radio source) and hot dust 
(MIR source) are intimately associated. We conclude that the 
emitting dust is located inside the H II region, mixed with the 
ionized gas. The A^-band spectra toward seven of the eight 
sources exhibit bright [Ne II] line emission, consistent with 
the association of these sources with ionized gas. 

The A^-band spectra show, for all eight observed sources, 
the presence of a broad and deep silicate absorption band 
at 9.7 /im, confirming their nature of deeply embedded ob- 
jects within massive dense cores of cold dust and molecu- 
lar gas. The H2 column densities derived from this absorp- 
tion feature are in the range (7- 17) x 10^^ cm"^, and are in 
good agreement (within a factor of ~ 2) with those estimated 
from 1.2 mm data (Paper II). Only toward one source, IRAS 
12383-6128, we detected the presence of PAH emission at 
8.6 and 11.3 /im. The lack of PAH emission in the spectra, 
with a small aperture and centered in the bright compact ob- 
ject, from the remaining sources is likely due to the presence 
of a hard radiation field which destroys these molecules in 
their neighborhood. 

We discovered five bright compact MIR sources which 
are not associated with radio continuum emission, and are 
thus prime candidates for hosting young massive protostars. 
In particular, objects IRAS 14593-5852 II (only detected at 
17.7 /im) and 17008-4040 I are likely to be genuine O-type 
protostellar objects. 



E.FE.M., D.M., G.G., K.J.B., and J.E.P gratefully ac- 
knowledge support from the Chilean Centro de Astroffsica 
FONDAP No. 15010003 and Proyecto BASAL PFB-06. 



REFERENCES 



Benjamin, R. A. et al. 2003, PASP, 115, 953 

Beuther, H., Churchwell, E. B., McKee, C. E, & Tan, J. C. 2007, in 

Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 165-180 
Beuther, H., Sridharan, T. K., & Saito, M. 2005, ApJ, 634, L185 
Bonnell, I. A. & Bate, M. R. 2002, MNRAS, 336, 659 
Bonnell, I. A., Bate, M. R., & Zinnecker, H. 1998, MNRAS, 298, 93 
Bonnell, I. A., Vine, S. G., & Bate, M. R. 2004, MNRAS, 349, 735 
Bronfman, L., Nyman, L.-A., & May, J. 1996, A&AS, 1 15, 81 
Brooks, K. J., Garay, G., Mardones, D., & Bronfman, L. 2003, ApJ, 594, 
L131 

Carey, S. J. et al. 2009, PASP, 121, 76 



Caswell, J. L. 1998, MNRAS, 297, 215 

Chapman, N. L., Mundy, L. G., Lai, S.-R, & Evans, N. J. 2009, ApJ, 690, 
496 

Cohen, M., Walker, R. G., Carter, B., Hammersley, P., Kidger, M., & 

Noguchi, K. 1999, AJ, 117, 1864 
De Buizer, J. M., Piiia, R. K., & Telesco, C. M. 2000, ApJS, 130, 437 
De Buizer, J. M., Radomski, J. T., Piiia, R. K., & Telesco, C. M. 2002, ApJ, 

580, 305 

De Buizer, J. M., Radomski, J. T, Telesco, C. M., & Piiia, R. K. 2005, ApJS, 
156, 179 

Draine, B. T. 2003, ARA&A, 41, 241 



A MULTIWAVELENGTH STUDY OF YOUNG MASSIVE STAR-FORMING REGIONS. III. 



15 



Engelbracht, C. W. et al. 2007, PASP, 1 19, 994 

Faundez, S., Bronfman, L., Garay, G., Chini, R., Nyman, L.-A., & May, J. 

2004, A&A, 426, 97 
Fazio, G. G. et al. 2004, ApJS, 154, 10 
Forster, J. R. & Caswell, J. L. 1989, A&A, 213, 339 

Garay, G., Brooks, K. J., Mardones, D., & Norris, R. R 2006, ApJ, 651, 914, 

(Paper 1) 

Garay, G. & Lizano, S. 1999, PASP, 111, 1049 

Garay, G., Mardones, D., Brooks, K. J., Videla, L., & Contreras, Y. 2007, 

ApJ, 666, 309, (Paper H) 
Home, K. 1986, PASP, 98, 609 

Jourdain de Muizon, M., Cox, R, & Lequeux, J. 1990, A&AS, 83, 337 
Kraemer, K. E. et al. 2003, ApJ, 588, 918 

Krumholz, M. R., Klein, R. I., & McKee, C. F. 2007, ApJ, 656, 959 
Mardones, D. 1998, Ph.D. Thesis, Harvard Univ. 
McKee, C. F & Tan, J. C. 2003, ApJ, 585, 850 

Meade, M. R. et al. 2005, GLIMPSE Legacy Science Data Products (ver 1.5; 

Madison: Dept. Astron., Univ. Wisconsin) 
Minier, V., EUingsen, S. R, Norris, R. R, & Booth, R. S. 2003, A&A, 403, 

1095 

Molinari, S., Pezzuto, S., Cesaroni, R., Brand, J., Faustini, F, & Testi, L. 

2008, A&A, 481,345 
Morris, P. W., Noriega-Crespo, A., Marleau, F. R., Teplitz, H. L, Uchida, 

K. I., & Armus, L. 2004, ApJS, 154, 339 
Natta, A. & Panagia, N. 1976, A&A, 50, 191 
Noriega-Crespo, A. et al. 2004, ApJS, 154, 352 
Ossenkopf, V. & Kenning, T. 1994, A&A, 291, 943 
Panagia, N. 1973, AJ, 78, 929 



Pascucci, I., Apai, D., Henning, T, Stecklum, B., & Brandl, B. 2004, A&A, 
426, 523 

Peeters, E. et al. 2002, A&A, 381, 571 
Povich, M. S. et al. 2007, ApJ, 660, 346 

Price, S. D., Egan, M. R, Carey, S. J., Mizuno, D. R., & Kuchar, T. A. 2001, 
AJ, 121, 2819 

Rathbome, J. M., Jackson, J. M., Chambers, E. T, Simon, R., Shipman, R., 

& Frieswijk, W. 2005, ApJ, 630, L181 
Reimann, H.-G., Linz, H., Wagner, R., Relke, H., Kaeufl, H. U., Dietzsch, 

E., Sperl, M., & Hron, J. 2000, in Proc. SPIE Vol. 4008, Optical and IR 

Telescope Instrumentation and Detectors, ed. M. lye & A. F. Moorwood, 

1132-1143 
Rieke, G. H. et al. 2004, ApJS, 154, 25 

Saviane, I. & Doublier, V. 2005, La Silla-Paranal Observatory. TIMMI2: 

User Manual (Issue 1.0; ESO) 
Siebenmorgen, R., Krugel, E., & Spoon, H. W. W. 2004, A&A, 414, 123 
Smith, H. A., Hora, J. L., Marengo, M., & Pipher, J. L. 2006, ApJ, 645, 1264 
Smith, M. D. & Rosen, A. 2005, MNRAS, 357, 1370 
Spitzer, L. 1978, Physical Processes in the Interstellar Medium (New York 

Wiley-Interscience, 1978. 333 p.) 
Walsh, A. J., Burton, M. G., Hyland, A. R., & Robinson, G. 1998, MNRAS, 

301, 640 

Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 
Werner, M. W. et al. 2004, ApJS, 154, 1 
Wood, D. O. S. & Churchwell, E. 1989, ApJ, 340, 265 
Zavagno, A., Cox, R, & Baluteau, J.-R 1992, A&A, 259, 241 



